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Abstract. We show that the radiation pressure of an intense optical field 
impinging on a perfectly reflecting vibrating mirror is able to entangle in a robust 
way the first two optical sideband modes. Under appropriate conditions, the 
generated entangled state is of EPR type [A. Einstein, et al., Phys. Rev. 47, 777 
(1935)]. 
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1. Introduction 

The radiation pressure acting on a movable mirror realizes an optomechanical 
coupling between the incident optical modes and the various vibrational modes of 
the mirror. The use of this coupling has been proposed many years ago for the 
implementation of quantum limited measurements of mechanical forces [1], as the 
interferometric detection of gravitational- waves [2], or atomic force microscopy [3]. 
Then optomechanical coupling has been proposed for quantum state engineering: for 
example, it has been shown that it may lead to nonclassical states of both the radiation 
field [4, 5], and the motional degree of freedom of the mirror [6]. 

The appearance of quantum effects in ponderomotive systems, paves the way 
for using them also for quantum information purposes [7]. In particular, quantum 
information can be encoded in the continuous quadratures of electromagnetic modes 
[8] , and many quantum communication protocols can be implemented using entangled 
states of optical fields [9]. A typical example of continuous variable entanglement 
is represented by the two-mode squeezed states at the output of a nondegenerate 
parametric amplifier [10]. However, it has been showed in Rcfs. [11, 12, 13] that 
radiation pressure can also be used to entangle two or more modes of a cavity with 
a movable and perfectly reflecting mirror. In the present paper, we get rid of the 
cavity and we show that radiation pressure provides a new way to entangle travelling 
electromagnetic modes. We shall consider the very simple case of an intense optical 
mode, incident on a single, vibrating, and perfectly reflecting mirror. A vibrational 
mode of the mirror will induce sideband modes on the reflected field, which will be 
shown to be entangled. Entanglement proves to be very robust with respect to the 
thermal noise acting on the mirror. In particular, for a well defined interaction time, 

1. e., a given time duration of the incident field, the two first sideband modes become an 
EPR-correlated pair, identical to the two-mode squeezed state, with the ratio between 
the optical and the mechanical frequency playing the role of the squeezing parameter. 

The paper is organized as follows. In Section II we shall derive the effective 
Hamiltonian of the system. In Section III we shall exactly solve the dynamics and 
provide a complete characterization of the reduced state of the system composed by 
the two reflected sideband modes. Section IV is for concluding remarks. 

2. Effective Hamiltonian of the system 

Let us consider a perfectly reflecting mirror and an intense laser beam impinging on 
its surface (see Fig. 1). For simplicity we consider only the motion and the elastic 
deformations of the mirror taking place along the spatial direction x, orthogonal to 
its reflecting surface. 

The electromagnetic field exerts a force on the mirror proportional to its intensity 
and, at the same time, it is phase-shifted by the mirror displacement from the 
equilibrium position [14, 15]. In the limit of small mirror displacements, and in the 
interaction picture with respect to the free Hamiltonian of the electromagnetic field 
and the mirror displacement field x(r,t) (r is the coordinate on the mirror surface), 
one has the following Hamiltonian [16] 




(1) 



where P(r, t) is the radiation pressure force [14]. All the continuum of electromagnetic 
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Figure 1. Schematic description of the system. A laser field at frequency loq 
impinges on the mirror oscillating at frequency Q. In the reflected field two 
sideband modes are excited at frequencies ui = u>o — Q and LU2 = <^o + fi- 



modes with positive longitudinal wave vector q and transverse wave vector k 
contributes to the radiation pressure force. Following Rcf. [14], and considering linearly 
polarized radiation with the electric field parallel to the mirror surface, we have 

Hr,t) = - A j dkj dqj dk' J dq'^L{M k ■ U k ,)u q 

x {a(k, q)a(k', q ) exp[-i(w + u')t + i(k + k') • r] 
+a f (k, g)a t (k', q 1 ) exp[i(w + - i(k + k') • r] 
+a(k, ^a^k', g') exp[-i(w - ui')t + i(k - k') ■ r] 
+^(k,q)a(k',q')aq>[i{Lj -u')t-i(k-k') ■ r}} , (2) 

where o(k, q) are continuous mode destruction operators obeying the commutation 
relations 

[a(k, q ),tf(k',q')] =6(k-k')S(q-q'). (3) 

Furthermore, the electromagnetic wave frequencies u> and u' are given by u> = 
yjc 2 (k 2 + q 2 ) and lo' = ^Jc 2 (k' 2 + q' 2 ) (c is the light speed in vacuum), and u,t, u q 
denote dimensionless unit vectors parallel to k, q respectively. 

The mirror displacement x(r, i) is generally given by a superposition of many 
acoustic modes [16]; however, a single vibrational mode description can be adopted 
whenever detection is limited to a frequency bandwidth including a single mechanical 
resonance. In particular, focused light beams are able to excite Gaussian acoustic 
modes, in which only a small portion of the mirror, localized at its center, vibrates. 
These modes have a small waist w, a large mechanical quality factor Q, a small effective 
mass M [16], and the simplest choice is to choose the fundamental Gaussian mode 



with frequency ft and annihilation operator b, with 



6, St 



1, 



' r(r - n 1 [ Se_<nt + Steint ] exp(-r 2 /™ 2 ). (4) 



H=- r » 2 
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By inserting Eqs. (2) and (4) in Eq. (1) and integrating over the variable r, one obtains 

fcr» V 2Wn I dk I dq / dk 7^^ (ufe ' Ufe ' H 

x {a(k, g)a(k / , <?') exp[-*(w + w')t - (k + k') 2 w 2 /4] 
+a t (k, g)a t (k', <?') exp[i(w + - (k + k') 2 w 2 /4] 
+a(k, ^(k', 9 ') exp[-i(w - J)t - (k - k') 2 w 2 /4] 
+a f (k, g)o(k / , cxp[i(w - w')t - (k - k') W/4] } 

x {be- int + Ve int } . (5) 

In common situations, the acoustical waist to is much larger than typical optical 
wavelengths [16], and therefore we can approximate exp { — (k ± k') 2 w 2 /4} w 2 /A-k ~ 
<5(k ± k') and then integrate Eq. (5) over k', obtaining 



// ■ — \/— — / </k / ,lq I dq'^L 



2n V 2Mfl 
x {— a(k, q)a(— k, q') exp[— i(uj + u)')t] 
— a^(k, q)af(— k, </) exp[i(w + u/)t] 
+a(k, g)a^(k, </) exp[— i(w — u/)t] 

+a^(k,q)a(k,q')cxp[i(LU-Lu')t}} x {fee"*"* + ftte*"*} . (6) 

We now make the Rotating Wave Approximation (RWA), that is, we neglect all the 
terms oscillating in time faster than the mechanical frequency f2. This means averaging 
the Hamiltonian over a time r such that fir 3> 1, yielding the following replacements 
in Eq. (6) 

2ir 

exp{±i<y ±w±ft)t} -> — 6{uj'±u±n). (7) 

T 

The parameter r is not arbitrary, but its inverse, 1/r = A^ et , is the detection 
bandwidth, that is, the spectral resolution of the detection apparatus employed. 

Since lu and lu' are positive and f2 is much smaller than typical optical frequencies, 
the two terms S(lu' + lu ± f2) give no contribution, while the other two terms can be 
rewritten as 

—6(u; , -u±n) = 2nAv det 5(q , -q ± )'^^-, (8) 



where q± = (lo ± fl) 2 /c 2 — k 2 . Integrating over q 1 we get 

J dk J dq-^j= |a(k, q)a) (k, q + ) b \Juj + Q 



H = - hAi> de 



2MVL 

+a(k, q)a) (k, <?_ ) SV w - fi 
+a t (k, g)a (k, q + ) tfy/uT+Q, 

+a t (k, <?)a (k, q_)Wu - ^} , (9) 

where we have used the fact that ou'(q±) — to ± f2. 

We now consider the situation where the radiation field incident on the mirror 
is characterized by an intense, quasi-monochromatic, laser field with transversal wave 
vector k , longitudinal wave vector q , cross-sectional area A, and power p. Since 
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this component is very intense, it can be treated as classical and one can approximate 
a(k, q) ~ a(k, g) in Eq. (9), where (with an appropriate choice of phases) 



a(k,g) = -i4/^|j(k-ko)%-?o), (10) 



with w = c^ko 2 + g 2 . 

Due to the Dirac delta, the only nonvanishing terms in the optomechanical 
interaction driven by the intense laser beam involve only two back-scattered waves, 
that is, the sidebands of the driving beam at frequencies wq ± fl, as described by 



H = iHAv det \ ——q 



P 



2MQ V hcj cA 



luq V w 



-a T (k , q+) b+J a T (k , ?_) 6 T 



w V 



u, " n -a(ko,g + )8t-J^a(k tol g_)SL > (ii) 



where now g± = ^/ (wo ± Sl) 2 /c 2 — fc 2 . The physical process described by this 
interaction Hamiltonian is very similar to a stimulated Brillouin scattering [17], even 
though in this case the Stokes and anti-Stokes component are back-scattered by the 
acoustic waves at reflection, and the optomechanical coupling is provided by the 
radiation pressure and not by the dielectric properties of the mirror. 

In practice, either the driving laser beam and the back-scattered modes are never 
monochromatic, but have a nonzero bandwidth. In general the bandwidth of the 
back-scattered modes is determined by the bandwidth of the driving laser beam and 
that of the acoustic mode. However, due to its high mechanical quality factor, the 
spectral width of the mechanical resonance is negligible (about 1 Hz) and, in practice, 
the bandwidth of the two sideband modes Av mo de coincides with that of the incident 
laser beam. It is then convenient to consider this nonzero bandwidth to redefine the 
bosonic operators of the Stokes and anti-Stokes modes to make them dimensionless, 



= ^ ^w, ;:k() (/ j = a^^acko.g.) (12) 



2irAv mode _ s Aq _ 

0,2 = V cA a q+ ' = V ~A a ^ 9+ ' ' 

so that Eq.(ll) reduces to an effective Hamiltonian 

H eff = -ihx(aib - a\&) - ih0(a 2 P - alb) , (14) 

where the couplings \ an d # are given by 



X = qvAv det 




2MO V Av mod JvjJ^ V ^0 



= C0S ^I 2M^A Vmode (15) 



X \IW^> (16) 
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with 0o = arccos(c(7o/w ), * s the angle of incidence of the driving beam. It is possible 
to verify that with the above definitions, the Stokes and anti-Stokes annihilation 



operators a\ and a 2 satisfy the usual commutation relations 



a, : , a 



3. System dynamics 

Eq. (14) contains two interaction terms: the first one, between modes di and b, 
is a parametric-type interaction leading to squeezing in phase space [10], and it is 
able to generate the EPR-like entangled state which has been used in the continuous 
variable teleportation experiment of Ref. [18]. The second interaction term, between 
modes d 2 and b, is a beam-splitter-type interaction [10], which may degrade the 
entanglement between modes di and b generated by the first term. In general, 
one has a system of three bosonic modes, coupled by a bilinear interaction. The 
corresponding Hamiltonian evolution of the system can be straightforwardly obtained. 
This Hamiltonian description satisfactorily reproduces the dynamics as long as the 
dissipative coupling of the mirror vibrational mode with its environment is negligible. 
This happens in the case of modes with a high-Q mechanical quality factor. In this 
case, the mechanical frequency f2 is sufficiently high (some MHz) so that the RWA of 
Eq. (7) can be made, and at the same time we can consider an interaction time, i.e., 
a time duration of the incident laser pulse, much smaller than the relaxation time of 
the vibrational mode (which can be of order of one second). 

The system dynamics can be easily studied through the (normally ordered) 
characteristic function $(/i, v, Q, where fi, v, ( are the complex variables corresponding 
to the operators a\,b, a 2 respectively. From the Hamiltonian (14) the dynamical 
equation for $ results 

i ( ***<9 d ± d 9 \ I 

$ = X [nv + ii *u* - n*— - li— -v*—- v— $ 
\ ov ov* Ofi a ii* J 

+ '(<*£ + <l --•£-'£)•• (17) 

with the initial condition 

$(i = 0) = cxp [-n\v\ 2 ] , (18) 

corresponding to the vacuum for modes di , a 2 and to a thermal state for mode b. The 
latter is characterized by an average number of excitations n = [exp^O/feeT) — 1] , 
T being the equilibrium temperature and ks the Boltzmann constant. Since the initial 
condition is a Gaussian state of the three-mode system, and the system is linear, the 
joint state of the whole system at time t is still Gaussian, with characteristic function 

*= ex P [-.%| 2 -£M 2 -£|CI 2 

+CV + Cii*v* + + ^*C* + ZX* + , (19) 

where 

A{t) = 2(r»-l)» { [1 " C ° S (26 * )] [ " (r2 _ 1} _ 1] 

+4r 2 [1 - cos(et)]} , (20) 

B{t) = / {[l + cos(2eO]rT(r 2 -l) + l-cos(2et)} , (21) 
Mr 1 ) 
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C W = 2(r2 _ 1)3 / 2 {^ 2sin ( e *) 

+ [n{r 2 - 1) - l] sin (26i)} , (22) 

V ^ = 2(r2 r i)3/ 2 { sin ( 9< ) + N r ' - !) " !] sin ( 20 *)} - (23) 

* (*) - 27^TF { [1 " cos (26 * )] [7T{r2 ~ 1} ~ l] 

+4 [1- cos (6*)]} , (24) 

^ = 2(/-l)2 { [1 ~ C ° S (29 * )] [ " (r2 ~ 1} ~ 1] 

+2(l + r 2 )[l-cos(9i)]} , (25) 

r = 0/x = [( w o + ^)/( w o — fi)] 1//2 and 6 = x(r 2 — l) 1 / 2 . The corresponding density 
operator can be expressed as 

P U .= I— ^^...CtJe-H'-M'-ICI' 

J TT J -K J -K 

XA(-/i)D 6 (-I/)^ 2 (-C), (26) 

where Z)j indicates the corresponding normally ordered displacement operator [19]. 

Here we are interested in the reduced state of the system composed by the two 
reflected optical sideband modes. They do not interact directly but their interaction 
is mediated by the optomechanical coupling of each mode with the moving mirror. 
This reduced state can be immediately obtained by tracing over the mirror mode b, 
obtaining 

$ 12 ( M , C) = cxp [-.4M 2 - £\(\ 2 + T (K + SO] ■ (27) 
Introducing the vector of field quadratures v = (Xi, Pi, X2, P2), where 
&j + d\ „ dj - a\ 

A >V' p ^^jt' j = 1 - 2 - <28) 

it is possible to connect the characteristic function (27) with the correlation matrix T 
of the two sideband modes, defined as Tij — (vjVj + VjVj}/2 — (vj)(vj), and which, 
in this case, is equal to 

( A + \ T \ 



r = 



A+\ -T 
T £ + \ 



(29) 



V -T £ + \ J 
The Gaussian state of the two optical modes is completely characterized by this 
correlation matrix. In particular we can employ the necessary and sufficient criterion 
for entanglement in the case of Gaussian states derived by Simon [20], which allows 
us to determine the parameter region where the two modes are entangled by the 
optomechanical interaction with the mirror vibrational mode. Simon's necessary 
and sufficient criterion for entanglement becomes, for the two-sideband modes state 
considered here, 

2 , ^2 



A 1\ 2 (£ 1 
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As it can be immediately seen from Eqs. (20)-(25), the dynamics of the system is 
determined only by three dimcnsionlcss parameters, Oi, r, and n, and it is periodic in 
<dt with period 2tt. The parameter r is determined by the ratio between the mechanical 
and the driving laser frequency, fl/wo and it is always very close to 1. For this reason 
in Fig. 2 we show the marker of entanglement T as a function of the scaled time Qt 
for < <dt < 2-7T, for different values of the initial mean thermal vibrational number 
n and at a fixed value of r, i.e., r — 1 = 2.5 x 10~ 7 . Fig. 2 clearly shows that the two 
sideband modes are almost always entangled (obviously except for t = 2kn, for integer 
k, due to the factorized initial condition) and that this entanglement is extremely 
robust with respect to the thermal noise acting on the mirror. In fact, the marker of 
entanglement T remains always negative, even at extremely large values of n. Fig. 2 
shows T for n — (full line), n = 10 5 (dashed line, almost coinciding with that for 
n = 0), n = 5 x 10 6 (dotted-dashed line), and n — 10 7 (dotted line). This result is 
extremely interesting because it shows that radiation pressure proves to be an efficient 
source of entangled travelling optical fields, not particularly affected by thermal noise. 
In this treatment we have only considered this kind of noise and we have neglected 
other more technical noise sources, such as the fluctuations of the intensity and the 
frequency of the incident driving field. In fact, the latter, differently from thermal 
noise, are negligible in today experiments involving optomechanical systems [21]. 

From Fig. 2, it is evident that T at exactly half period, <dt — it, is unaffected by 
thermal noise, i.e., it is independent of n. In fact, using Eqs (20)-(25) and (30), it is 
easy to see that T(Qt = n) = — 4(r + r 3 ) 2 / (r 2 — l) 4 (actually T has been scaled exactly 
by the modulus of this value in Fig. 2). This independence of n is a manifestation 
of quantum interference between the dynamical effects of the two interaction terms 
in Eq. (14). It is therefore interesting to see the reduced state of the two reflected 
sideband modes just at this interaction time, <dt = tt. Using Eqs. (20)-(25), one has 
the following correlation matrix 

r(et = 7r)= (3i) 
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1 

2 




2 



2r(l+r 2 ) 








2r(l+r 2 ) 
(r 2 -l) 2 





2r(l+r 2 ) 

( r 2_!)2 



4r 2 , 1 
(r 2 -l) 2 + 2 



V 







4r 2 , 1 
(r 2 -!) 2 + 2 / 



It is easy to check that this is just the correlation matrix of an EPR-like entangled 
state, i.e., identical to the two-mode squeezed state generated by a parametric amplifier 
[10]. The correspondence between the squeezing parameter £ and our parameters is 
sinh£ <-> 2r(l + r 2 )/(r 2 — l) 2 , so that, in practice one can reach very high two-mode 
squeezing by choosing r — 1 — > 0, i.e., by decreasing the ratio CI/uj . This fact can 
be equivalently seen by evaluating the variances of the linear combination of field 
quadratures, typical of the two-mode squeezed state, i.e., 



The time behavior of the variances (32) and (33) is depicted in Fig. 3. It is worth 
noticing that the typical EPR correlations are available for a wide time range around 
Qt = ir. Furthermore, they are robust against thermal noise since the curves are 
shown for n = 10 5 , and they practically coincide with those for n = 0. That is in 
agreement with the results shown in Fig. 2. 

4. Results and Conclusions 

In conclusion, we have seen that, under appropriate conditions, the radiation pressure 
of an intense optical pulse impinging on a perfectly reflecting mirror couples efficiently 
the first two sideband modes induced by a vibrational mode of the mirror. This 
optomechanical interaction actually entangles the two modes, and quite surprisingly, 
the resulting entanglement is extremely robust with respect to the thermal noise acting 
on the mirror. In particular, due to quantum interference, if the interaction time, 
i.e., the duration of the incident driving pulse, is appropriately tailored, the reduced 
state of the reflected sideband modes is a two-mode squeezed state. Therefore, rather 
unexpectedly, radiation pressure could be employed as a new source of entangled two- 
mode squeezed state, in which the role of the squeezing parameter is played by the 
ratio between the optical and the mechanical frequency u>q /CI. These interesting results 
hold if the time duration of the pulse satisfies two conditions: i) it has to be much 
longer than the period of mechanical oscillations, because we have to average over 
these oscillations; ii) it has to be shorter than the vibrational mode relaxation time 
(mechanical damping has been neglected in our treatment). These conditions are 
not extremely easy but they could be realized using micro-opto-electro-mechanical- 
systems (MOEMS) for example [22]. Possible values are: damping rates 7 m ~ 1 Hz, 



p = 10W, luq ~ 2 x 10 15 Hz, CI ~ 5 x 10 8 Hz, Ais det ~ 10 7 Hz, Av mode ~ t' 1 ~ 10 3 
Hz, and M ~ 10~ 10 Kg, yielding x - $ - 5 x 10 5 Hz, and 6 ~ 10 3 Hz. Therefore, 



A_ = ((X, - X 2 f) = ((P, + P 2 f) , 
A + = {{X l +X 2 f) = {{P l -P 2 f). 



(32) 
(33) 



For Qt — ir, one has 




(34) 



(35) 
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Figure 3. The quantities A+(©i) (a) and A— (©<) (b) are plotted vs the scaled 
time ©t. The values of other parameters are r = 1 + 2.5 X 1CT 7 and n = 10 5 . 

using well tailored and long pulses (order of 1 ms), one could realize this new source 
of two mode squeezing. 
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